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Abstract 

Chinese honeybee Apis cerana (Ac) is one of the major Asian honeybee species for local apiculture. However, Ac is frequently 
damaged by Chinese sacbrood virus (CSBV), whereas Apis mellifera [Am) is usually resistant to it. Heterospecific royal jelly 
(RJ) breeding in two honeybee species may result in morphological and genetic modification. Nevertheless, knowledge on 
the resistant mechanism of ;Am to this deadly disease is still unknown. In the present study, heterospecific RJ breeding was 
conducted to determine the effects of food change on the larval mortality after CSBV infection at early larval stage. 2-DE and 
MALDI-TOF/TOF IVIS proteomic technology was employed to unravel the molecular event of the bees under heterospecific 
RJ breeding and CSBV challenge. The change of Ac larval food from RJC to RJM could enhance the bee resistance to CSBV. 
The mortality rate of Ac larvae after CSBV infection was much higher when the larvae were fed with RJC compared with the 
larvae fed with RJIVl. There were 101 proteins with altered expressions after heterospecific RJ breeding and viral infection. In 
Ac larvae, 6 differential expression proteins were identified from heterospecific RJ breeding only, 21 differential expression 
proteins from CSBV challenge only and 7 differential expression proteins from heterospecific RJ breeding plus CSBV 
challenge. In Am larvae, 17 differential expression proteins were identified from heterospecific RJ breeding only, 26 
differential expression proteins from CSBV challenge only and 24 differential expression proteins from heterospecific RJ 
breeding plus CSBV challenge. The RJM may protect Ac larvae from CSBV infection, probably by activating the genes in 
energy metabolism pathways, antioxidation and ubiquitin-proteasome system. The present results, for the first time, 
comprehensively descript the molecular events of the viral infection of ;Ac and Am after heterospecific RJ breeding and are 
potentially useful for establishing CSBV resistant populations of Ac for apiculture. 
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Introduction 

Apis mellifera {Am) and A. cerana (Ac) are major honey bee 
species in the global beekeeping industry [1,2]. They are heavily 
infected by different vital viruses [3,4]. Chinese sacbrood virus 
(CSBV) is the most stricken pathogen of Ac, which results in severe 
and deadly infections of the colony and eventually losses of the 
entire colony [5]. As a small RNA virus (picorna-like virus) that 
has an icosahedral virion with a diameter of 26-30 nm, CSBV 
genome consists of a single positive-strand RNA molecule with 8.8 
kilo bases (kb) [7]. Although CSBV shows a close genetic 
relationship to its western counterpart sacbrood virus (SBV), there 
is no cross infection [8]. Since this viral disease broke out in 1972 
in southern China [8,9], some eflForts have been made to study this 
virus, such as diagnostic methods (electron microscopy, enzyme- 
linked immunosorbent assay and reverse transcription-polymerase 
chain reaction (RT-PCR)) [10], and the control of this disease by 
RNA interference [11]. It is interesting that A. mellifera is not 
sensitive to CSBV in general beekeeping practice [12]. Similar to 
all other insects, the honeybees lack a classically adaptive immune 



system as in the case of mammalian [6] . To survive, they have 
evolved cellular and humoral immune responses to cope with 
microbial infections. A. mellifera may develop cellular and 
humoral immune responses to various pathogens such as bacteria 
[13,14], viruses [15], microsporidian [16-18] and Varroa mites 
[19,20]. Usually, the viral infection causes cell apoptosis, tissue 
damage and even functional disorder of the whole organism and 
all these changes can be further reflected in proteome alteration 
[21]. Recent study revealed the pathological mechanism of 
Chinese sacbrood disease to Ac [6]. However, no information 
on the difference in the immune responses of two bee species to 
CSBV is reported. 

In the honey bees, queen and workers have different behavior 
and reproductive capacity despite possessing the same genome. 
The primary substance that leads to this differentiation is royal 
jelly (RJ), which contains a range of proteins, carbohydrates, lipids, 
minerals, vitamins, and a large number of bioactive substances, 
especially immunological peptides and antibacterial proteins [22— 
24]. Major Royal Jelly Proteins (MRJPs) are the prime RJ 
ingredients, which are crucial in regulating reproductive matura- 
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tion [22]. There are quantitative differences in nucleic acids and 
protein composition in fresh RJ between Am and Ac [24] . The 
most recent discovery is that RJ contains microRNAs which may 
play a role in caste differentiation [25-27]. Whether the BJ also 
influences the resistant behavior oiAm to CSBV is unknown. The 
molecular event for CSBV resistance of Am needs to be also 
determined. 

Two-dimensional gel electrophoresis (2-DE) based proteomics 
technology and matrix assisted laser desorption ionization time-of- 
flight mass spectrometry (MALDI-TOF MS) are powerful 
molecular tools. In recent years, proteomic techniques have been 
successfully used to profile the proteome change in insect growth 
and development [6,28-30]. The present study employed gel- 
based (2-DE) and shotgun proteomic (label-free LC-MS based) 
strategies, which have complementary natures, to gain an in-depth 
understanding of the resistant mechanism of Am to the fatal CSBV 
disease by comparison of the proteome-wide change of the healthy 
and sick worker larvae of both and^c fed with RJ either from 
Am (RJ^/[) or Ac (RJC). 

In the present study, we performed comparison analyses of the 
proteome in Ac and Am after heterospecific RJ breeding (i.e. Am 
larvae fed with RJC or Ac larvae fed with RJM), and then CSBV 
challenge, by two-dimensional electrophoresis and MALDI-TOF 
MS analysis, to systematically search for the molecular basis of the 
crossbreeding and CSBV challenge. A total of 101 differentially 
expressed nonredundant proteins (>3 folds change) were identi- 
fied. Furthermore, gene expression patterns in two bee species 
were investigated at mRNA levels. The results show that 
differential expression proteins are involved in the regulation of 
metabolism and response to CSBV challenge. 

Materials and Methods 

Ethics Statement 

Honey bee colonies (Am and Ac) were raised at Conghua, 
Guangdong Province, China (1 13" 17' E. 23°8" N), by Guangdong 
Entomological Institute, by standard beekeeping techniques. 
According to the present regulation from Ministry of Agriculture 
Am and Ac used in this studies are not endangered or protected 
species. No specific permissions were required for performing 
these experiments. 

Honey Bees 

To obtain age controlled second instar larvae, the queen was 
caged on a comb and left to lay eggs for 6 h. Twenty hours after 
larval eclosion, or 92 h after oviposition, the comb containing 
second instar larvae was retrieved from the colony, and placed in 
the laboratory for treatments at 32-34°C. All the larvae used in 
this study were detected by Reverse Transcription-Polymerase 
Chain Reaction (RT-PCR) method for the absence of the 
following viruses, black queen cell virus (BQCV), chronic bee 
paralysis virus (CBPV), deformed wing virus (DWV), kashmir bee 
virus (KB\^, Chinese sacbrood virus (CSB\'^ and Israeli acute 
paralysis virus (lAPV), with the primers from Ai et al. [4]. RNA 
was extracted using Trizol Reagent (Invitrogen, Carlsbad, USA) 
and cDNA synthesis was performed using PrimeScript 1st Strand 
cDNA Synthesis Kit (TaKaRa, Kyoto, Japan). Amplification 
profile of PGR consisted of an initial 2-min denaturation at 94°G, 
followed by 35 cycles of 30 s at 95°G, 30 s at 55°C, 1 min at 72°C 
and finally 7-min cycle at 72"C. No signs of clinical American 
foulbrood or other viral disease were observed in these larvae. 



Virus 

GSBV-infected larvae with typical symptom were collected from 
an apiary at Xinhua Village from Hunan Province, and k(-pt at - 
80°C for less than 10 days befor(; use. The presence of CSBV virus 
was confirmed by observation of the morphological symptom 
under electronic microscopy and RT-PCR according to Chen et 
al. and Yan et al. [3,10]. 

To obtain CSBV, the larvae (approximately 0.2 g) infected by 
CSBV were ground in 6 mL sterile phosphate-buflFered solution 
(PBS) (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HP04, 1.8 mM 
KH2PO4, pH 7.4) with a sterile grinder. The resulting solution 
was centrifuged at 14 000 g at 4°C for 10 min, and the 
supernatant was further passed through a 0.45 ^m cell filter first, 
then through a 0.22 |J.m cell filter [31]. The CSBV concentration 
was quantitated by absolute quantification assay [11]. CSBV 
polyprotein gene (SBVl) was amplified to detect copy number of 
CSBV according to Liu etal. [11]. PGR reactions were carried out 
in triplicate in an Mx3000P Real-Time PGR System (Stratagene, 
California, USA), using SYBR_Gr(-cn fBrilliant II SYBR-Green 
QPGR Master Mix; Stratagene, California, USA). The CSBV 
concentration for the experiments was approximately 
7.9x10^ copies/mL. 

Harvest of RJM and RJC 

BJM and RJC were produced according to standard practices 
in China [32]. Briefly, the queen was confined inside a queen 
excluding cage. Queen cups with young larvae (one-day old) were 
introduced into the colony and allowed to be fed by workers for 2 
days. 3-day old larvae were first carefully removed by using either 
a grafting tool or a pair of forceps, then the royal jelly was removed 
by using a spatula. 

Larvae Fed with RJM or RJC and Challenged with CSBV 

To determine the effects of the changing food (RJs) on the 
mortality and protein expression of the third instar larvae, which 
are the most sensitive stage to CSBV infection, the second instar 
larvae of Ac and Am were collected carefully from the combs, 
respectively transferred to 96-cell culture plates containing 30 (iL 
RJC or RJM, in each cell, then placed in a cabinet (Sanyo, Tokyo, 
Japan) at 80% relative humidity at 32-34°C. After 24 h, the tiiird 
instar larvae of both bee species fed with RJM or RJC were 
transferred to new plates, and challenged with 30 |J,L CSBV 
solution per cell, prepared as described above. After 8 h, the larvae 
were washed with sterile distilled water for 3 times, and kept at — 
80°C for use. 8 treatments (Table 1) were established with 3 
replicates containing 20 larvae in each replicate. Each 20 larvae 
were pooled as one biological replicate and three independent 
biological replicates were produced. 

The same treatments were also estabUshed to examine the 
mortality of bee larvae fed with RJM or RJC and challenged by 
CSBV after 24 h. The death of the larvae was confirmed when 
they had no response after stimulating with a soft tip, and by the 
color change of their bodies (from white to yellow or even dark). 

Protein Extraction and 2-DE 

Larval protein extractions were carried out according to the 
previously described method with some modifications [33,34]. 

Bee larvae were manually homogenized for 30 min and 
sonicated for 2 min on ice in lysis buffer (8 M urea, 0.2% w/v 
Bio-Lyte 3/10 Ampholyte (Bio-Rad, USA), 4% CHAPS, 65 mM 
DTT) containing Protease Inhibitor Cocktail (Calbiochem, 
Germany)(w/1 = 1 mg/4 |iL). The mixture was centrifuged 15 
000 g at 4°C for 30 min, and the supernatant was collected. The 
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Table 1. Treatments of the two honey bee species by RJs and CSBV. 



No. 


Honey bees 


Treats 


Abbreviation 


A 


A. cerana 


fed on A. cerana royal jelly 


AC-RJC 


B 




fed on A. cerana royal jelly, and challenged with CSBV 


AC-RJC+CSBV 


C 




fed on A. mellifera royal jelly 


AC-RJM 


D 




fed on A. mellifera royal jelly, and challenged with CSBV 


AC-RJM+CSBV 


E 


A. mellifera 


fed on A. cerana royal jelly 


AlVl-RJC 


F 




fed on A. cerana royal jelly, and challenged with CSBV 


AM-RJC+CSBV 


G 




fed on A. mellifera royal jelly 


AM-Rm 


H 




fed on A. mellifera royal jelly, and challenged with CSBV 


AM-RJIVl+CSBV 



doi:1 0.1 371 /journal.pone.Ol 02663.t001 



resulting precipitates were treated three times with the same 
procedures. All the supeniatants were centrifuged 15 000 g at 4°C 
for 60 min. The final supernatant was used or stored at — 80°C for 
2-DE. Protein concentrations were determined by the Bradford 
method using Modified Bradford Protein Assay Kit (Sangon, 
Shanghai, China). 

The 2-DE was performed according to the methods described 
previouslyand the manufacturer's instruction [35]. The first 
dimension (isoelectric focusing) was conducted using the IPGphor 
lEF system (Bio-Rad, California, USA) at 20°C. For analytical 
gels, 350 |J,g protein was solubHized in 400 |J.L rehydration 
solution (8 M urea, 0.2% w/v Bio-Lyte 3/10 Ampholyte, 4% 
CHAPS, 65 mM DTT, 0.001% w/v bromophenol blue), and 
loaded onto a 17 cm pH 3-10 NL IPG strip (Bio-Rad). Focusing 
was performed for 13 h at 50 V, 1 h at 500 V, 1 h at 1000 V, and 
30 min and 6 h at 8000 V (total = 48 kVh). The IPG strips were 
equilibrated as previously described. The second dimension was 
performed with 13% (w/v) SDS-polyacrylamide gels using the 
Protean II xi 2D MulticeU system (Bio-Rad, California, USA). 
Proteins were stained with silver nitrate, and gels were digitized 
using Image Scannerll (Amersham Biosciences, USA). Digitized 
2-DE gel patterns were edited and matched using the PDQUEST 
V8.0.1 software package (PDI, Humingtoii Station). Triplicate 
experiments were run to confirm the reproducibility of results. 

Protein Identification by MALDI-TOF/TOF MS 

Protein digestion and peptide extraction were done according to 
our previously established protocol [36] . Spots of interest in gels 
staining with silver nitrate were cut out, washed, reduced, S- 
alkylated with iodoacetamide and in-gel digested at 37°C 
overnight with sequencing grade porcine trypsin (Promega, 
USA). After extraction in extractant of 50% ACN (Fisher, 
Waltham, MA,USA) and 2.5% TFA (Sigma, St Louis, MO, 
USA), peptide mixtures were analyzed using a saturated solution 
of 5 mg/ mL a-cyano-4-hydroxycinnamic acid (CHC A, Sigma, St 
Louis, MO, USA) in ACN containing 0. 1 % TFA (Sigma, St Louis, 
MO, USA) (50/50 v/v) using a 4800 Proteomics Analyzer 
equipped with matrix assisted laser desorption ionization time- 
of-flight mass spectrometry (MALDI-TOF MS) (Applied Biosys- 
tems, Framingham, MA, USA). For MS calibration, the trypsin 
autolysis peptides were used as internal calibrants. Monoisotopic 
peak masses were automatically determined within the mass range 
of 800^000 Da, with a minimum S/N of 50. Five of the most 
intense ion signals were selected as precursors for MS/MS 
acquisition. Combined MS and MS/MS queries were performed 
with the MASCOT search engine (V2.1, Matrix Science, UK) 
embedded in GPS-Explorer Software (V3.6, Applied Biosystems, 



Framingham, MA, USA), using the A. mellifera database (Gene 
DB). MASCOT protein scores (based on combined MS and MS/ 
MS spectra) of greater than 61 were considered statistically 
significant (P^0.05). The individual MS/MS spectrum with 
statistically significant (confidence interval >95%) best ion score 
(based on MS/MS spectra) were also accepted. 

Based on the UniProt Knowledgebase (http://www.uniprot. 
org/), identified proteins were submitted to KEGG (Kyoto 
Encyclopedia of Genes and Genomes) and classified according 
to KEGG pathway maps (http:/ /www.genome.jp/kegg/pathway. 
html) and Gene Ontology (http://www.geneontology.org/), 
grouped on the basis of their biological process of GO terms 
(http://www.blast2go.com) [19]. Protein interaction networks of 
the differentially regulated proteins were analyzed using the online 
database resource Search Tool for the Retrieval of Interacting 
Genes (STRING: http://string-db.org) [37,38]. 

Quantitative Real-Time PCR (qRT-PCR) 

To verify the gene expression of differentially expressed 
proteins, 11 proteins (FBA, LOC408516, MDH, SOD, Jafracl, 
LOC725646, AGOl, Che-3, Aid, Pxd and Tpi) were randomly 
selected for qRT-PCR analysis of their encoding genes by 
quantitative real time PCR, according to our previously method 
with some modifications [19]. Briefly, total RNA was isolated 
using Trizol reagent (Invitrogen, Carlsbad, USA). The RNA was 
dissolved in diethylpyrocarbonate (DEPC) treated water and 
quantified by measuring ultraviolet absorbance at 260 and 
280 nm. Reverse transcription was performed on 1 [ig of total 
RNA using PrimeScript RT reagent Kit with gDNA Eraser 
(Takara, Kyoto, Japan). The gene-specific primers for the analysis 
are listed in Table 2. Expression oi aclin gene (GI: 406122) was 
used as an internal control. The efficiency of each primer set was 
first validated by constructing a standard curve through five serial 
dilutions. PCR reactions were carried out in triplicate in an 
Mx3000P Real-Time PCR System (Stratagene, California, USA), 
using SYBRGreen (SYBR Premix Ex Taq II, Takara, Kyoto, 
Japan). A control without template was included in all batches. 
The PCR program began with a single cycle at 95°C for 10 min, 
40 cycles at 95°C for 15 s and 48°C for 60 s. Afterwards, the PCR 
products were heated to 95°C for 15 s, cooled to 48°C for 15 s 
and heated to 95 °C for 15 s, in order to measure the dissociation 
curves and to determine a unique PCR product for each gene. 
mRNA levels were calculated relative to actin expression using the 
Mx3000P Software (version 4.1) (Agilent, California, USA). The 
fold change was calculated using the 2 '^^'''^ method. Each sample 
was analyzed independently and processed in triplicate. 
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Table 2. Primers for qRT-CR analysis of 12 differentially expressed genes. 



Product n3iTi6 


Gene ID 


Forw3rd priiner seciuence 


Reverse primer sec|uence 


Amplicon size (bp) 


Actin 


gi|406122 


ATGCCAACACTGTCCTTTCTGG 


GACCCACCAATCCATACGGA 


150 


FBA 


gi|66526635 


TTTCTAGTATCTAAAGCATG 


ACTGCTATTGCTACTGCT 


228 


LOC408516 


gi|66524882 


GATCCAGATGCACCAAG 


GGACAGAATCGTGGAAA 


235 


MDH 


gi|66513092 


AAGGCTGGCACAGGTTC 


TAAATGCAGCGATCCCA 


239 


SOD 


gipjuoy 1 U4 


AL 1 1 Li 1 1 1 LLu 1 u 1 A 


/- Ar~ A 1 1 \f~ A A^~r'^~ ATTA 
LALA 1 1 1 LAALLLA 1 1 A 


229 


Jafraci 


gi|66548188 


AAACTCATTGCAGCATC 


GGGAGGTCTGTTGATGA 


250 


LOC725646 


gi|1 10755974 


ACCCGAATTGTACTTTA 


TAGATGAAGGTGCTCAA 


221 


AGOl 


gi| 110777044 


TGGCCCAGATCAAGTAGAGC 


AATTTGATAGCGTTTGTGGTGAT 


200 


Che-3 


gi|l 10757336 


TCCTGTCCGAATTTTTACCTGT 


GAAGCTGCGTTTGCGTCTA 


250 


Aid 


gi| 110748949 


TGCGTACTGTTCCACCT 


GAGGGCTAAGGCTAACA 


233 


Pxd 


gi| 110757934 


CGCTGTTCCAAGAGGCT 


ACGACAAACGCCAGATC 


250 


Tpi 


gi| 148224276 


TCAGTAACAGCAGGAAAT 


TTTATGTCATGTTCCTGCAA 


269 


doi:1 0.1 371 /journal.pone.Ol 02663.1002 



Western Blot 

To further verify the variation tendency of differentially 
expressed proteins identified by the proteomic approaches, 
glyceradehyde- 3 -phosphate dehydrogenase 2 isofome l(GAPDH), 
fmarylacetoacetase hydrolase (FAH), heat shock protein (Hsp)60, 
phosphoglycerate Idnase (PGK), triosephosphate isomerase 1 
(Tpi), CuZn superoxide dismutase (SOD) and fructose-bispho- 
phate aldolase (FBA) were randomly selected for Western Blot 
analysis by the method of Han et al. [6] with minor modifications. 
Briefly, equal amount of protein sample (30 |tg/lane) were 
separated by stacking (5%) and separating (12%) SDS-PAGE 
(sodium dodecyl sulfate polyacrylamide gel electrophoresis) gels 
and then transferred to a nitrocellulose transfer membrane 
(0.2 |lm pore size) (Invitrogen, Eugene, OR) using the iBlot 
apparatus (Invitrogen, Eugene, OR). After blocking, the mem- 
branes were incubated over night at 4°C with primary rabbit 
polyclonal antibodies of anti-GAPDH, FAH, Hsp60, PGK, Tpi, 
SOD, and FBA antibodies (Abeam, Cambridge, MA) at a dilution 
of 1:1000. Following three washes, the membranes were further 
incubated with horseradish peroxidase-conjugated rabbit anti-goat 
secondary antibody (BOSTER, Wuhan, China) at a dilution of 
1:2000 for 1 h. Immunoreactive protein bands were detected using 
the DAB Western Blotting Substrate (BOSTER, Wuhan, China) 
and quantified by densitometry using Quantity-one image analysis 
system (Bio-Rad, Hercules, CA). Actin was detected simultaneous- 
ly as loading control of the analysis, which immunoreactived by 
BCIP/NBT Western Blotting Substrate (BOSTER, Wuhan, 
China). Quantification of the protein bands was conducted by 
scanning the films and importing the images into the Quantity 
One 1-D analysis software (Bio-Rad, Hercules, CA). Scanning 
densitometry was used for semiquantitative analysis of the data. 

Data Analysis 

The statistical analysis of gene expression was performed by 
one-way ANOVA test followed by Tukey's test for pair 
comparisons (SPSS version 16.0, SPSS, Inc.). Comparison of 
larval mortality differences between different groups was per- 
formed by independent samples t test. An error probabiUty p< 
0.05 was considered statistically significant. 



Results 

Effects of RJM and RJC on the Mortality of Two Bee 
Species after CSBV Challenge 

Experiments were carried out to examine whether RJM would 
enhance the resistant ability o{ Ac larvae to CSBV infection. The 
larval mortalities of the treatments after 56 h were indicated in 
Figure 1. The dead larvae showed the symptom of CSBV 
infection. Without CSBV infection, the mortalities of Ac larvae 
fed with RJC or RJM were not significantly different (t = 2.040, 
df= 4, jb = 0.1 1 1; Figure 1 A); meanwhile, the mortalities oi Am 
larvae fed with RJM or RJC were also not significantly different 
(t= 1.572, df=4, j5) = 0.191; Figure 1 E). This indicated that the 
RJs did not significantly influence the mortalities of the larvae of 
two bee species, at least at key early larval stage. 

However, the mortality (72.2%) of larvae fed with RJC and 
subsequent CSBV infection was significantly higher than that of 
the larvae fed with RJC (26.7%) but without CSBV (t= -9.642, 
df= 3.801, p = 0.001; Figure 1 B); When Ac larvae were fed with 
RJM and subsequently challenged by CSBV, their mortality was 
not significantly different than that of those larvae fed with RJM 
and without CSBV (t= -2.447, df=4, = 0.071; Figure 1 C); 
Interestingly, when Ac larvae were fed with BJC or RJM and 
subsequently challenged with CSBV, the mortality (72.2%) of 
larvae fed with RJC was significantiy higher than that of those 
larvae fed with RJM (27.8%,) (t= 10.674, df=4, j!) = 0.000; 
Figure 1 D). This indicated that RJM conferred antiviral activity 
for Ac larvae. 

When Am larvae were fed with RJC or RJM, and subsequendy 
challenged with CSBV, no significant differences in mortalities of 
Am larvae were observed (t= 1.572, df= 4, jf) = 0.191, Figure 1 E; 
t=-0.941, df=4, /) = 0.400, Figure 1 F; t=-0.512, df=4, 
^ = 0.636, Figure 1 G; t = 2.000, df = 4, = 0.1 16, Figure 1 H). 
So the BJs did not influence the mortality of Am larvae after 
CSBV challenge. As usual, Am is not sensitive to CSBV infection. 

Qualitative Comparisons and Identification of 
Differentially Expressed Proteins 

To determine the effects of the RJ on the protein expression of 
two bee species after CSBV infection, the protein spots were 
visualized by 2-DE from the bee larvae of Ac and Am fed with 
RJM or RJC, and with RJM or RJC for 24 h and then challenged 
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Figure 1 . The larval mortalities of A. cerana and A. mellifera under different treatments after 56 hours. RJC: royal jelly from A. cerana; FUM: 
royal jelly from A. mellifera; AC-RJC: 3-day Ac larvae fed with RJC. AC-RJM: 3-day Ac larvae fed with RJM. AC-RJC+CSBV; 3-day Ac larvae fed with RJC, 
and infected with CSBV. AC-RJM+CSBV: 3-day Ac larvae fed with RJM, and infected with CSBV. AM-RJC: 3-day Am larvae fed with RJC. AM-RJM: 3-day 
Am larvae fed with RJM. AM-RJC+CSBV: 3-day Am larvae fed with RJC, and infected with CSBV. AM-RJM+CSBV: 3-day Am larvae fed with RJM, and 
infected with CSBV. The data are averages from three replicates. The error bars indicate standard deviations. Bars with different letters mean 
significantly different (P<0.05, independent samples t test). 
doi:1 0.1 371/journal.pone.01 02663.g001 



by CSBV for 8 h. 8 h after CSBV challenged, there was no dead 
larvae, although the virus was detected in the challenged larvae. 

The total proteins were separated on 2-DE gels spanning pH 3- 
10, silvcT stained, and analyzed by MS. Protein levels were 
expressed as percentage volume, which corresponds to the 
percentage ratio between the volume of a single spot and the 
total volume of aU spots present in a gel. The mean values of spot 
intensity were calculated using at least three gels. Spots showing 
more than 15% variation were not considered (Student's test, with 
7 degrees of freedom, P<0.()5). Little deviation was observed in 
the patterns on replica gels. 

530-670 spots were revealed on the silver-stained 2-DE patterns 
of the total proteins from dilferent larvae, depending on bee 
species and the treatments (Figure 2). Protein spots were 
distributed over the 3-10 pH range. 

Quantitatively, 296 total spots showed significant changes of 
expression (by a factor of s3-folds and j!><0.05). AU dilferentially 
expressed spots were identified by a MAIDI-TOF/TOF mass 
spectrometer, and submitted to the NCBI nr database or the 



NCBI EST_ database. 163 differentially expressed proteins were 
successfully identified (Figure 2), among which 44 were identified 
as MRJPs, and 1 19 were clustered to 101 non-redundant proteins. 

The up-regulated and down-regulated proteins from the 3-day old 
Ac and^m larvae fed with RJC or RJM, and with CSBV infection 
were listed (Table 3). 

In Ac larvae, 6 differential expression proteins (6 spots) were 
identified from heterospecific RJ breeding only, 21 differential 
expression proteins (22 spots) from CSBV challenge only and 7 
differential expression proteins (7 spots) from heterospecific RJ 
breeding plus CSBV challenge. In Am larvae, 17 differential 
expression proteins (17 spots) were identified from heterospecific 
RJ breeding only, 26 differential expression proteins (36 spots) 
from CSBV challenge only and 24 differential expression proteins 
(31 spots) from heterospecific BJ breeding plus CSBV challenge 
(Figure 3). 

Some proteins were not identified either because their 
abundance was too low to produce enough spectra or because 
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Figure 2. 2-DE pictures of A. cerana and A. mellifera 3-day old larvae fed with RJs and by CSBV infection. RJC: royal jelly from A. cerana; 
RJM: royal jelly from A. mellifera; A: 3-day old Ac larvae fed with RJC. B: 3-day old Ac larvae fed with RJM. C: 3-day old Ac larvae fed with RJC and by 
CSBV infection. D: 3-day old Ac larvae fed with RJIVI and CSBV infection. E: 3-day old Am larvae fed with RJC. F: 3-day old Am larvae fed with RJIVl. G: 3- 
day old Am larvae fed with RJC and CSBV infection. H: 3-day old Am larvae fed with RJM and by CSBV infection. 
doi:l 0.1 371/journal.pone.01 02663.g002 



the database search scores can not yield unambiguous results (> 

95%). 

GO (Gene Ontology) Functional Term Enrichment 

To understand the functions of difiFerentiaUy expressed proteins 
from the bee larvae, which resulted from heterospecific RJ 
breeding or CSBV infection, 101 identified proteins were analysed 
by gene ontology (GO) mapping respertively. At least one GO 
term could be assigned to 88 of 101 detected difiFerentiaUy 
expressed proteins (Table 3). Thirteen (12.8%) proteins were 
deemed of unknown function after blast searches against the nr 
databases and annotation augmentation. The GO terms assigned 
to individual proteins were shown in Table 3. Based on the GO 
annotations, proteins were grouped by their broad biological 
processes (BP) (Figure 4, 5). 

Three major BPs of the up-regulated proteins from the Ac 
larvae fed with RJM were metabolic processes (40.0%), cellular 
processes (20%) and biological regulation (20%), whereas those of 
down-regulated proteins were mainly metabolic processes (17.6%), 
cellular processes (17.6%) and developmental process (11.8%). 
The major BPs of the increasing proteins from Ac fed with BJG, 
and then by CSBV challenge were metabolic process (18.2%), 
cellular processes (18.2%) and response to stimulus (18.2%), 
whereas the major BPs of the decreasing proteins were metabolic 
processes (37.5%), cellular processes (33.3%) and biological 
regulation (12.5%). The major BPs of the increasing proteins 

40 n 

35 - 

.1 30 - 




from Ac fed with RJM, and then by CSBV challenge were 
metabolic processes (25.0%), cellular processes (25.0%) and 
cellular component organization or biogenesis (16.7%), whereas 
the major BPs of the decreasing proteins were metabolic processes 
(25.0%), cellular processes (25.0%) and biological regulation 
(12.5%). 

Three major BPs of the increasing proteins from the Am larvae 
fed with RJM were metabolic processes (22.2%), cellular processes 
(11.1%) and cell killing (11.1%), whereas those of decreasing 
proteins were metabolic processes (32.0%), cellular processes 
(28.0%) and response to stimulus (16%). The major BPs of the 
increasing proteins from Am fed with RJM, and then by CSBV 
challenge were cellular processes (28.6%), biological regulation 
(19.0%) and cellular component organization (9.6%), whereas the 
major BPs of the decreasing proteins were metabolic processes 
(37.5%), cellular processes (29.2%) and response to stimulus 
(16.7%). The major BPs of the increasing proteins from Am fed 
with RJC, and th(-n by CSBV challenge were metabolic processes 
(60.9%), cellular processes (iOA"/o) and response to stimulus 
(4.3%), whereas the major BPs of the decreasing proteins were 
locahzation (28.5%), response to stimulus (21.4%) and metabolic 
processes (21.4%). 

Network Analysis of Differentially Expressed Proteins 

Proteins function jointly in a hving cell through networks by 
forming protein- protein interactions (PPI), modifications and 




Figure 3. Numbers of up-regulated and down-regulated proteins from the bee larvae fed with RJs and challenged by CSBV. On the 

basis of triplicate replications analyses, only the proteins that changed >3.0-fold in relative ratios (p<0.05) were considered. RJC: royal jelly from A. 
cerana; RJM: royal jelly from A. mellifera; AC-RJM: 3-day old Ac larvae fed with RJM. AC-RJC+CSBV: 3-day old Ac larvae fed with RJC and infected with 
CSBV. AC-RJM+CSBV: 3-day old Ac larvae fed with RJM and infected with CSBV. AM-RJC: 3-day old Am larvae fed with RJC. AM-RJC-fCSBV: 3-day old 
Am larvae fed with RJC and infected with CSBV. AM-RJM-fCSBV: 3-day old Am larvae fed with RJM and infected with CSBV. 
doi:1 0.1 371/journal.pone.01 02663.g003 



PLOS ONE I www.plosone.org 



12 



August 2014 I Volume 9 | Issue 8 | el 02663 



Heterospecific Feeding and CSBV Challenging 



GO:006S007 
biological regulation. 
20K 



COCX}OS152 
metabolic process,. 
40K 




GO:0071840 
cellular component 
organization or 
biogenesis, 20% 



60«X)9987 
cellular process, 20K 



GO0050S96 
response to 
stimulus, 11.80%. 

GO:006S007 
biological 
gulation, llSWiV 

GO:0032SO: 
multicellular 
;anismal process.^ 
11.80K 

GO«0OSi52 
metabolic process. 
17 60* 




G0«I71840 
cellular component 

organization or 
biogenesis. 11.80% 

00:0032502 



developmental 
process. 11.80% 



GO:0009987 cellular 
process, 17.«0K 



G0:0040011 
ocomotion. 5.80% 



GO:00S1179 
localization. 4.5% 
GO:000S2S3 
cell proliferation. 4 5% 

GO:0000003 
reproduction. 4.5%. 

GO 0065007 biological 
regutation,9.1% 

GO:0071840 cellular, 
component 
organization or 
biogenesis, 4.5% 



GO 0032501 
multicellular 
organismal process 




response to 

stimulus, 18.2% 



GO:0032502 
developmental process. 
9.1% 



GO0O08152 
metabolic process, 
18.2% 



GO:00099S7 
cellular process 18.2% 



B 



GO:006S007 
biological 
regulation, 8.3% 

GO.00325O1 
multicellular 
organismal process, 
83% 



GO:0008152 
metabolic process. 



GO 0051179 
localization, 8.3% 



GO.0071840 
cellular component 
organization or 
biogenesis, 16.7% 




GO:0009987 
el lular process; 



GO.t)0S089« 
response to 
stimulus, 8 30%. 

00:0051179 
localization, 4.20%. 



00 0065007 
biological 
gulation, 12.50%. 



GO:0008152 
metabolic process, 
37.SOK 



GO4)0S089« 
response to stimulus, 
1430% 

00:0001906 
cellkilling. 14 30%. 

00:0032501 
multicellular 
organismal proces: 
14.30% 




00:0071840 
cellular component 
organization or 
biogenesis, 4.20% 



00:0009987 
cellularprocess, 
33 30% 




00:0032502 
developiYwntal 
process. 14.30% 



000009987 
ellular process, 
14.30% 



_0O«)51704 
multi-organism 
process, 14.30% 



00 0065007 
biological regulation. 
14.30% 



Figure 4. Based on the GO annotations, the proteins from A. cerana larvae fed with RJM or RJC, and by CSBV challenge, were 
grouped by the broad biological processes (BP). A, B, C: up-regulated proteins; D, E, F: down-regulated proteins, RJC: royal jelly from A. cerana; 
RJM: royal jelly from A. mellifera; A, D: 3-day /Ac larvae fed with RJIVl; B, E: 3-day /Ac larvae fed with RJC, and infected with CSBV; C, F: 3-day ,Ac larvae fed 
with RJIVl, and infected with CSBV, 
doi:10,1371/journal.pone.0102663,g004 



regulation of expression relationships. Of all the identified proteins 
in this study, 107 proteins with the estabUshed functions based on 
Drosophila melanogaster were selected and their PPI mapping was 
constructed using STRING. The result showed that all the treats 
almost divided into two distinct PPI clusters with several others 
operating without connectivity. The first cluster primarily consist- 
ed of cytoskeletal proteins, while the second cluster mainly 
consisted of proteins involved in metabolism (Figure 6). 

Verification of Differentially Expressed Proteins 

To test the tendency of protein expression between its encoding 
gene at the transcript level, 1 1 proteins (FBA, LOC408516, MDH, 
SOD, Jafracl, LOC725646, AGOl, Che-3, Aid, Pxd and Tpi) 
were randomly selected for qRT-PCR analysis of their encoding 
genes. The trend of mRNA expression of the selected genes 
showed consistent patterns with their protein expression (Figure 7). 



Western Blot Analysis 

Western blot analysis was conducted to verify the expression of 
the proteins which are closely related to energy metabolism of 
heterospecific royal jelly and CSBV resistance, such as GAPDH, 
PGK, Tpi, FAH, SOD, and FBA. The results are consistent well 
with the proteomic data (Figure 8). 

Discussion 

Heterospecific RJ breeding means during the growth or 
reproduction period of the colony or individual, changing the 
comb from Ac to Am or from Am to Ac or feeding RJ artificially. In 
recent years, heterospecific RJ breeding between Ac and Am has 
been concerned in the bee industry [39] . This technology was used 
between different species of western honeybees in 1957. Am fed 
with RJC showed some morphological characteristic of^c, such as 
proboscis, forewings, abdomen and thorax, and body weight of 



PLOS ONE I www.plosone.org 



13 



August 2014 | Volume 9 | Issue 8 | el 02663 



Heterospecific Feeding and CSBV Challenging 




Figure 5. Based on the GO annotations, the proteins from A. 
mellifera larvae fed with RJM or RJC, and by CSBV challenge, 
were grouped by the broad biological processes (BP). A, B, C: 

up-regulated proteins; D, E, F: down-regulated proteins. RJC: royal jelly 
from A. cerana; RJM: royal jelly from A. mellifera; A, D: 3-day Am larvae 
fed with RJC; B, E: 3-day Am larvae fed with RJM, and infected with 
CSBV; C, F: 3-day Am larvae fed with RJC, and infected with CSBV. 
doi:1 0.1 371 /journal.pone.01 02663.g005 

new adults [23,40], and even influenced the genotype of malate 
dehydrogenase II [41]. It was reported that Am could be resistant 
to varroa mites when fed with RJC [42-44]. Interestingly, the 
present results showed that the change of larval food from RJC 
to RJM could enhance Ac resistance to CSBV, under the 
laboratory condition. The mortality of Ac larvae after CSBV 
infection was much lower when they were fed with RJM compared 
with those fed with RJC. However, the change of RJM to RJC did 
not lower the viral resistance of Am larvae. Whether BJM 
stimulates the viral resistance of Ac by inducing the antiviral 
protein expression or RJM molecules can directly inhibit the viral 
replication in bee larvae needs further study. 

To determine whether the proteins induced by the food change 
are probably involved in the bee resistance to CSBV, 2-DE and 
MALDI-TOF-MS based proteomic strategies were used. Accord- 
ingly, 101 nonredundant proteins were identified as being 
differentially expressed among Ac and Am larvae by heterospecific 
RJ breeding and CSBV challenge. 

Differentially Expressed Proteins from Hererospecific RJ 
Breeding 

The results showed that protein expression was influenced by 
the RJ. Within 6 differentially expressed proteins from 3-day old 
Ac larvae fed with RJM, 3 proteins CG6987-PA (SF2), vasa 
intronic gene (Vig), and GAPDH were up-regulated, and 3 
proteins (Act5C, ben, Che-3) down-regulated. SF2 is a RNA 
recognition motif, a highly abundant domain in eukaryotes found 
in proteins involved in post-transcriptional gene expression 
processes [45]. Vig has been identified as a component of 
Drosophila RISC involved in RNA interference [46] . GAPDH is a 
key enzyme in the glycolytic pathway. So the up-regulated proteins 
are associated with gene expression process, RNA interference and 



glycolytic metabolism. Of the down-regulated proteins, Gamma- 
actin (Act5C) is an ubiquitous protein involved in the formation of 
filaments that are a major component of the cytoskeleton. 
Ubiquitin conjugating enzyme E2 (E2) is an ubiquitin conjugating 
enzyme [47,48] . The ubiquitin proteasome system enables cells to 
selectively recognize and degrade damaged or potentially harmful 
proteins [49]. Abnormal CHEmotaxis family member (Che-3) is 
one of the AAA-H superfamily. In Caenorhabditis elegans Che-3 is 
specifically responsible for the retrograde transport of the 
anterograde motor, kinesin-II, and its cargo within sensory in 
chemosensory neurons [50]. Therefore, the down-regulated 
proteins are related with the cytoskeleton arrangement, protein 
metabolism and cellular processes. 

Correspondingly, 17 differentially expressed proteins from 3- 
day old larvae of Ain fed with RJC included 6 up-regulated 
(pORF2, yeUow-e2, MDH, Oat, RpS12, one unknown protein), 
and 11 down-regulated proteins (Hsc70-5, Tpx-1, Ef2b, FAH, 
invA, Cct5, Rpsl2, Hsp60,Gapdhl, HexllO, one unknown 
protein). For the up-regulated proteins, pORF2 is a retrovirus 
reverse transcriptase with RNA-directed DNA polymerase activity 
[39]. Yellow-e2 responses to fung-us/G"^/G~, and functions in 
caste determination, and is influenced by environmental factors 
[51]. MDH is the key enzymes in the citric acid cycle, and possess 
L-malate dehydrogenase activity elevated by oxidative stress [52]. 
Ornithine aminotransferase precursor (Oat) possesses ornithine- 
oxo-acid transaminase activity in metabolic process [53]. For the 
down-regirlated proteins, T-complex Chaperonin 5 (Cct5) is 
involved in productive folding of proteins [54] . The methylation 
of the HexllO gene in A. mellifera is regulated at the 
developmental stage and in a caste-dependent manner [55]. 
HSP60 and Hsc70-5 are heat shock proteins responding to stress. 
Generally, HSPs act as molecular chaperones faciUitating protein 
maturation in cells and activating innate and acquired immune 
responses [6]. They are excellent immunomodulators against a 
wide variety of pathogens in protecting host [56]. Thioredoxin 
peroxidase 1 (Jafracl) is in response to oxidative stress [57]. 
Elongation factor 2 (Ef2b) is a GTPase, involved in the 
translocation of the peptidyl-tRNA [58]. Fumarylacetoacetase 
hydrolysis (FAH) is responsible for the hydrolysis of fumarylace- 
toacetate to generate fumarate and acetoacetic acid [59] . Putative 
invasion protein (InvA) is a prominent inner-membrane compo- 
nent of the Salmonella type III secretion system (T3SS) apparatus, 
which is responsible for regulating virulence protein export in 
pathogenic bacteria [60]. 

From the above description, differentially expressed proteins 
from heterospecific RJ breeding were related to metabohsm 
mainly in TCA and glycolytic pathway and response to microbe 
and oxidative stress. RJ protein difference between /4c and Am was 
reported [24]. The majority of the identified proteins were major 
royal jelly proteins (MRJPs) with MRJPl being the most 
abundant. Peroxiredoxin 2540, glutathione S-transferase SI, and 
MRJP5 were detected only in the RJ of ^4. mellifera ligustica, and 
MRJP7 was found only in the RJ of^. cerana cerana [24]. Maybe 
different components of two royal jellies changed the metabolism 
pathways of honey bees, conferring to the resistant difference to 
viral infection. 

Differentially Expressed Proteins from Ac Larvae Fed with 
RJC and by CSBV Challenge 

Comparing to the 3-day olAAc larvae fed with RJC and without 
CSBV infection, 21 altered proteins were detected from 3-day old 
Ac larvae fed with RJC and CSBV infection, in which 9 proteins 
being up-regulated (Vha55,Hsp70, RpS3A, MTAl-like, pxd, 
AGOl, phosphatidylethanolamine-binding protein, HsplO, an 
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PFPs 

RpS4 0.999 
RpS20 0999 
RpL27A 0.999 
RpSSa 0999 
RpL8 0 999 
KpSS 0 999 
RpLH 0 999 
KpS13 0999 
RpS24 0 999 
sop 0999 



B 




PFPs 

Oap<fli2 0 999 
Pjjt 0 999 
OipdhI 0 999 
Eno 0999 
RpSlA 0 999 
Pgfc 0 999 
C09961 0 999 
Mhc 0 999 
RpLZI 0 998 
PyK 0998 




PFTi: 

RpS3 0.999 
RpS13 0 999 
RpLll 0 999 
Rf>S23 0.999 
RpSlS 0999 
RpL30 0 999 
RpLlOAb 0999 
RpLZ3 0 999 
RpS2A 0 999 
RpLlBA 0999 



Figure 6. Protein-protein interaction networlcs of identified differentially expressed proteins from nutritional crossbreeding and 
CSBV challenge were constructed by Search Tool for the Retrieval of Interacting Genes (STRING). Lines represent the existence of the 
different types of evidence used in predicting the associations. A red line indicates the presence of fusion evidence; a green line, neighborhood 
evidence; a blue line, co-occurrence evidence; a purple line, experimental evidence; a yellow/ line, text mining evidence; a light blue line, database 
evidence; a black line, coexpression evidence. PFPs: Predicted Functional Partners. A: 3-day Ac larvae fed with RJM. B: 3-day Ac larvae fed with RJC, 
and infected with CSBV. C: 3-day Ac larvae fed with RJM, and infected with CSBV. D: 3-day Am larvae fed with RJC. E: 3-day Am larvae fed with RJM, 
and infected with CSBV; F: 3-day Am larvae fed with RJC, and infected with CSBV. 
doi:10.1371/journal.pone.0102663.g006 
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Figure 7. Verification of 11 differentially expressed proteins at 
mRNA level by qRT-PCR analysis. The normalized relative mRNA 
levels (>1) of the genes of the corresponding differentially expressed 
proteins indicate up-regulation and those (<1) indicate down- 
regulation. Error bar is standard deviation. The names of the proteins 
are referred in Table 2. 
doi:1 0.1 371 /journal.pone.01 02663.g007 

unknown protein), and 12 proteins down-regulated 
(V12B01_09746, Aats-tyr, Aid, AKR-1, RpS9, Jafracl, PGK, 
Hsp60, mid, wph, LOC552007, LOC409705). 

The up-regulated proteins were mainly involved in the stress 
response and in the defense from viral infection. Such as, vacuolar 
h (Vha55) couples ATP hydrolysis to the build up of a H+ 
gradient. Metastasis associated gene 1 (MTAl) is part of the 
NURD (nucleosome remodeling and deacetylating) complex and 
plays a role in cellular transformation and metastasis [61]. 
Peroxidase-like (Pxd) is response to oxidative stress. Argonaute 1 
(AGOl) is part of the RNA-induced silencing complex (RISC) in 



the RNA interference pathway which can defense insects from 
virus infection [62,63]. 

The down-regulated proteins were mainly involved in the 
energy process and in the development regulation. For example, 
tyrosyl-trna cytoplasmic (Aats-tyr) are involved in amino acid 
synthesis, cell cycle control, RNA shear modification and transport 
of tRNA [64,65] . Fructose-bisphosphate aldolase (FBA) and aldose 
reductase-like isoform 1 (AKR-1) are involved in glycolysis. PGK 
and pyruvate kinase are involved in energy metablism. PGK is the 
enzyme responsible for the first ATP generating step of glycolysis 
and the in vivo activation of 1-nucleoside pro-drugs effective 
against retroviruses such as HIV and hepatitis [66]. Pyruvate 
kinase (PK) regulates the final rate-limiting step of glycolysis. In 
tumor cells, the glycolytic pyruvate kinase isoenzyme M2 (PKM2, 
M2-PK) determines whether glucose is converted to lactate for 
regeneration of energy or used for the synthesis of cell buUding 
blocks [67]. Cell proliferation only proceeds when metabolism is 
capable of providing a budget of metabolic intermediates that is 
adequate to ensure both energy regeneration and the synthesis of 
cell building blocks in sufficient amounts. In this study, the three 
important enzymes of glycolysis were significantly down regulated 
after Ac larvae fed with BJC and by CSBV challenge, indicating 
that the energy metabolism pathway may be important in antiviral 
activity of the bees. Prohibitin protein wph (Wph) is required for 
larval metabolism or for the progression of the larva into a pupa. 
LOC409705 is responsible for imaginal disc-derived wing 
morphogenesis. These proteins were significantly down-regulated, 
probably direcdy impacting the larval development. 

As indicated in Figure 1 , the Ac larvae fed with RJC were not 
resistant to CSBV infection. Although the up-and down-regulated 
protein expressions were detected from those bee larvae, it seemed 
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Figure 8. Western blot Analysis of differentially expressed proteins. Six proteins were selected. Proteins were separated on a 1 2% SDS-PAGE 
gel and immunoblotted with anti-GAPDH, anti-PKG2, anti-Tpi, anti-FAH, anti-SOD, and anti-FBA. p-actin was used as an internal control. (A) The 
Western blot images of GAPDH, PKG2, TPI, FAH, SOD, FBA and Actin. The protein spots from Figure 2. (B) Quantification of target protein 
expression.The relative fold change of GAPDH, PKG2, Tpi, FAH, SOD, FBA (normalized by Actin) was visulized by using Quantity One Software. The 
expression of GAPDH, SOD, FBA and Tpi were up-regulatedd. On the other hand, FAH and PGK were down-regulated. Error bar is standard deviation. 
doi:1 0.1 371/journal.pone.01 02663.g008 
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that these proteins were not actively involved in effectively 
protecting the larvae fi-om CSBV infection. Whether some Ac 
expressed proteins were active for helping viral infection needs 
further study. 

Differentially Expressed Proteins from Ac Larvae Fed with 
RJM and by CSBV Challenge 

Comparing to the 3-day old^c larvae fed with RJM but without 
CSBV infection, 4 highly expressed protein spots (cIF-5A. Tpi. 
CG2330RA, Act88F) and 3 down-regulated proteins (CG10703, 
HSP60, yellow-e2) were found from the Ac larvae fed with RJM 
and by CSBV infection. 

Eukaryotic translation initiation factor 5a (eIF-5A) plays an 
important role in protein translation extending, and in stress 
tolerance [68]. Triosephosphate isomerase (Tpi) is a glycolytic 
enzyme [69] . Act88F (actin, indirect flight muscle-like) is involved 
in phagocytosis [51]. HSP 60 and yeUow-e2 are both activated by 
environmental factors as discussed above. They were down- 
regulated may indicated that the influence of CSBV were 
weakened when the Ac larvae fed with RJM. 

Interestingly, CSBV infection did not significantly increase the 
mortality of larvae fed with RJM comparing to those fed with 
RJM and without viral challenge (Figure. 1). So RJM may protect 
Ac larvae from viral infection probably by promoting energy 
metabolism and activating phagocytosis. Another possibility is that 
the specific components from RJM, such as Peroxiredoxin 2540, 
glutathione S-transferase SI, and MRJP5, which are detected only 
in the RJ of A. mellifera ligustica [24], may inhibit the viral 
replication. 

Differentially Expressed Proteins from Am Larvae Fed 
with RJC and by CSBV Challenge 

In comparison to 3-day old^?re larvae fed with RJC and without 
CSBV infection, 24 proteins were detected from 3-day old Am 
larvae fed with EJC and CSBV infection, in which 17 proteins 
were up-regulated (Cryll (lambda-crystallin homolog), 
Dper\GL19229, Cat, CG9117 (metallo-beta-lactamase), 
LOC725919 (transposase), Eif, RpS20, serine protease easter 
(Sp7), LOC551766 (atp synthase subunit mitochondrial), SOD, 
Abnormal wing disks protein (Awd), FBA, ornithine aminotrans- 
ferase (Oat), GADPH, larval cuticle protein a3a-like (CPR31A)) 
and 7 proteins down-regulated (Hsp6(), Beta-Spec, Stm-Mlck, 
Yellow-h, FABP, NPC2, an unknown protein). 

Most of the up-regulated proteins were involved in energy 
metabolism, such as Cryll, CG9117, LOC725919, Sp7, 
LOC551766, FBA, Oat, LOC410122, and GAPDH. The 
regulation pattern of FBA and GADPH was opposite to that in 
the Ac larvae fed with RJC and chaflenged with CSBV, maybe 
due to different species of honey bees. But in response to viral 
invasion, two bee species showed the same way in regulation of 
antioxidative enzymes, such as SOD, CuZn superoxide dismutase 
and Catalase, which were biomarkers to evaluate the toxic effects 
of organisms [70]. The proteins (Awd and CPR31A) related to 
development were up-regulated, to some degree, but the 
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cytoskeletal proteins (Beta-Spec and Stm-Mlck) were significantly 
down-regulated. These differentiafly expressed proteins might be 
involved in protecting Am larvae from CSBV infection, although 
these larvae were fed with RJC (Figure 1). Am larvae are not 
sensitive to CSBV infection [6]. The food change from RJM to 
RJC did not lower the resistant ability of Am to CSBV. 

Differentially Expressed Proteins from Am Larvae Fed 
with RJM and by CSBV Challenge 

In response to CSBV chaUenge, 11 up-regulated proteins 
(SS1G_11931, Pdi, FBA, CG5590, Hsc70-3, PGK, Cp309, 
Hsc70-4, GAPDH, Rpn9, Mhc), and 15 down-regulated proteins 
(RpS3, Jafracl, CG6084, H.sp60, Aats-tyr, ALD, Hsp70, Lsp2, 
Vha68, LOC552007, HexllO, NPC2, Roel, LOC725646) were 
found from the larvae fed with RJM and with CSBV infection. 

Among the up-regulated proteins, those involved in energy 
metabolism were the majority, including PDI, FBA, ALD, PGK 
and GAPDH. This was in according with the expression pattern of 
Ac larvae fed with RJM. 26S proteasome (Rpn9) plays a 
fundamental role in eukaryotic homeostasis [71,72]. The ubiqui- 
tin-proteasome system mediated viral protein degradation consti- 
tutes a host defense process against some RNA viral infections 
[73] . In the present study Rpn9 and HSPs were both up-regulated, 
indicating the possible involvement of these proteins in antiviral 
activity. Correspondingly, chaperonins or store proteins were 
down-regulated. This may guarantee the normal growth of Am 
larvae. As a source of amino acids for tissue reconstruction during 
pupal development, hexamerins are first synthesized by a larval fat 
body and released into the hemolymph where they accumulate to 
extraordinarily high concentrations to bufld adult structures [74] . 

Conclusions 

Although Ac is frequentiy damaged by CSBV, Am is usually 
resistant to this virus. The present study showed that the change of 
Ac larval food from RJC to RJM could enhance the bee resistance 
to CSBV, at least at early larval stage. Proteomic technology was 
employed to unravel the molecular events of the bees under 
heterospecific RJ breeding and CSBV chaUenge. There were 101 
proteins with altered expressions after heterospecific RJ breeding 
and viral infection in two bee species. The RJM may protect Ac 
larvae from CSBV infection, probably by activating the genes in 
energy metabolism pathways, antioxidation and ubiquitin-protea- 
some system. The results, for the first time, comprehensively 
descript the molecular events of the viral infection of and Am, 
and are potentially useful for establishing CSBV resistant 
populations of Ac for apiculture. 
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